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Abstract 

This  paper  reports  experimental  and  numerical  investigations  on  the  combustion  characteristics  of  a 
lean  methane-air  mixture  in  a  heated  porous  sand  bed.  The  porous  bed  consisted  of  sand  (Si02)  particles 
with  a  mean  particle  diameter  of  0.56  mm.  The  horizontally  placed  quartz  tube  was  heated  externally  to 
initiate  the  combustion  reaction  in  the  porous  bed  combustor.  The  stabilized  flame  location  curve  as  a 
function  of  averaged  mixture  velocity  was  obtained  for  various  equivalence  ratios.  Contrary  to  the  earlier 
finding  of  a  C-shape  flame  stabilization  behavior,  a  new  S-shape  behavior  was  observed  in  the  present 
study.  This  can  be  divided  into  three  regimes:  high,  moderate,  and  low  velocity  regimes.  In  the  low  velocity 
regime,  flame  with  very  weak  luminosity  was  confirmed  and  the  stabilized  flame  location  moved  down¬ 
stream  with  the  increase  of  the  mixture  velocity.  For  the  moderate  velocity  regime,  a  stable  flat  flame 
was  observed  and  the  flame  location  moved  upstream  with  the  increase  of  the  mixture  velocity.  An  oscil¬ 
latory  flame  behavior  was  observed  in  the  high  velocity  regime.  In  this  oscillatory  mode,  the  flame  front 
oscillated  with  a  characteristic  time  period  of  the  order  on  1  h  and  increased  with  the  increase  of  the  mix¬ 
ture  velocity.  In  order  to  further  understand  these  experimental  results,  one-dimensional  computational 
studies  with  detailed  chemistry  and  heat  transfer  mechanisms  were  carried  out.  The  computational  results 
were  in  good  agreement  with  experimental  observations.  The  computations  showed  that  solid-to-solid  radi¬ 
ation  played  a  significant  role  in  the  flame  stabilized  location.  From  the  examination  of  the  flame  structure, 
it  was  found  that  the  flame  behavior  in  the  low  velocity  regime  was  similar  to  that  of  the  flameless  com¬ 
bustion  mode. 
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1.  Introduction 

Fluidized  bed  combustion  systems  have  been 
utilized  for  industrial  incinerators  because  of  their 
inherent  advantage  of  low  NO*  emissions  and 
higher  combustion  efficiency  [1,2].  These  systems 
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use  waste  and  low  calorific  value  fuels  such  as  sew¬ 
age  sludge  or  livestock  waste  as  primary  fuel. 
Since  such  waste  contains  large  amounts  of  water 
and  incombustible  material,  it  is  difficult  to  sus¬ 
tain  stable  combustion.  Due  to  this  reason,  gas¬ 
eous  fuels  such  as  methane  or  propane  are 
introduced  to  enhance  the  combustion  stability 
in  a  fluidized  bed.  The  usage  of  an  external  electric 
heater  is  also  proposed  to  improve  the  combus¬ 
tion  stability  in  such  incinerators.  The  combustion 
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Nomenclature 

A 

Cross-sectional  area  of  the  combus- 

CQk 

Molar  rate  of  production  of  kt h 

tion  tube 

species 

At 

Outer  surface  area  of  tube/unit  length 

s 

Porosity 

cv 

Specific  heat  of  gas 

p 

Density 

Dr 

Particle  diameter  of  spheres 

<p 

Equivalence  ratio 

h 

n  gS 

Volumetric  heat  transfer  coefficient 

G 

Stefan-Boltzmann  constant 

between  gas  and  solid 

kr 

Radiative  thermal  conductivity 

Nondimensional 

2 

Thermal  conductivity 

Nu 

Nusselt  number 

% 

Radiative  heat  flux  in  axial  direction 

Pr 

Prandtl  number 

T 

Temperature 

Re 

Reynolds  number 

Th 

Heater  temperature 

u 

Mixture  velocity 

Subscripts 

Yk 

Mass  fraction  of  kt h  species 

g 

Gas  phase 

Vk 

Diffusional  velocity  of  Mi  species 

s 

Solid  phase 

Wk 

Molecular  weight  of  kt h  species 

k 

Species  index 

phenomenon  in  a  fluidized  bed  combustor  is 
intrinsically  correlated  with  both  the  gas  reaction 
in  bubbles  and  the  gas  reaction  in  the  condensed 
phase  [3].  Although,  the  gas  reaction  in  bubbles 
can  be  treated  based  on  typical  gas  combustion 
phenomena,  the  gas  reaction  in  the  condensed 
phase  of  the  fluidized  bed  is  remarkably  different 
from  that  of  ordinary  gas  combustion.  To  under¬ 
stand  the  gas  reaction  in  the  condensed  phase, 
porous  media  is  an  appropriate  choice  due  to  its 
simplicity  compared  to  the  complicated  combus¬ 
tion  process  in  fluidized  bed  combustion  systems. 

Combustion  in  porous  media  is  also  known  as 
filtration  gas  combustion  and  has  been  investigat¬ 
ed  by  many  researchers  [4-10]  to  understand  the 
combustion  characteristics  in  the  condensed 
phase.  Katani  et  al.  [11]  have  studied  the  stability 
limits  and  combustion  characteristics  in  a  bundled 
ceramic  tube  burner.  They  observed  that  the  flame 
is  located  in  the  downstream  or  upstream  region 
or  within  the  combustion  tube  depending  on  the 
flow  rates  and  equivalence  ratios.  Similar  combus¬ 
tion  regimes  have  been  observed  by  Babkin  et  al. 
[12-14].  They  observed  two  different  velocity 
regimes  for  the  stabilized  flame  location  in  a  por¬ 
ous  bed:  the  low  velocity  regime  (LVR)  and  the 
high  velocity  regime  (HVR).  In  the  low  velocity 
regime,  the  reaction  zone  and  temperature  profiles 
in  gas  and  solid  phases  propagate  as  a  united 
structure  with  a  velocity  on  the  order  on  lcm/s. 
In  the  high  velocity  regime,  the  gas  combustion 
zone  moves  almost  independently  of  the  solid 
temperature  profile.  In  between  these  two  regimes, 
transient  flame  propagation  was  observed.  How¬ 
ever,  the  behavior  of  stabilized  flame  within  the 
porous  media  was  not  fully  investigated. 

Various  researchers  [15-17]  have  emphasized 
the  importance  of  various  heat  transfer 


mechanism,  especially  the  role  of  solid-to-solid 
radiation.  Chen  et  al.  [15]  have  reported  that  the 
solid-to-solid  radiation  used  to  preheat  an 
unburned  mixture  leads  to  an  increase  of  burning 
velocity  and  reduced  flame  thickness.  Further 
study  of  the  importance  of  solid-to-solid  radiation 
has  been  conducted  by  Yoshizawa  et  al.  [16].  They 
simulated  combustion  in  porous  media  by  using 
three  different  sections  of  the  burner  and  conclud¬ 
ed  that  solid-to-solid  radiation  plays  an  important 
role  in  excess  enthalpy  burning.  In  spite  of  the  fact 
that  the  earlier  works  investigated  the  effect  of  sol¬ 
id-to-solid  radiation,  the  effect  of  solid-to-solid 
radiation  on  the  stabilized  flame  location  remains 
unclear. 

In  the  present  study,  the  combustion  character¬ 
istics  of  a  fuel-lean  methane-air  mixture  in  a  heated 
porous  sand  bed  were  investigated  both  experimen¬ 
tally  and  numerically.  The  different  flame  stabiliza¬ 
tion  regimes  were  classified  into  appropriate 
categories  and  the  flame  behaviors  in  these  regimes 
were  discussed.  The  effect  of  solid  radiation  on  the 
flame  stabilization  was  also  examined  by  using  dif¬ 
ferent  models  with  and  without  solid-to-solid  radi¬ 
ation  heat  transfer.  The  flame  behavior  was 
examined  based  on  analysis  of  the  flame  structure 
in  each  corresponding  regime. 

2.  Experimental  setup 

A  schematic  diagram  of  the  experimental 
apparatus  is  shown  in  Fig.  1.  The  sand  was 
packed  in  a  quartz  tube  as  a  porous  material. 
The  thermal  conductivity  and  the  permeability 
(for  the  single  phase  flow)  of  the  sand  bed  were 
2.7  x  10  1  W/mK  and  2.0  x  10  7  cm2,  respectively 
[18].  The  mean  diameter  of  Si02  sand  particle  was 
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Fig.  1.  Experimental  setup. 


0.56  mm,  which  was  determined  by  sieve  analysis. 
The  combustion  tube  had  an  inner  diameter  of 
27  mm  and  a  length  of  450  mm.  An  external  elec¬ 
tric  heater  was  used  to  heat  the  porous  sand  bed. 
The  heater  was  located  1 50  mm  downstream  of 
the  point  where  the  mixture  entered  the  porous 
sand  bed.  The  initial  temperature  distribution  of 
the  porous  sand  bed  is  shown  in  Fig.  2.  The  high¬ 
est  temperature  in  the  porous  sand  bed  was  main¬ 
tained  at  around  1330  K. 

Lean  methane-air  mixtures  with  equivalence 
ratios  of  0.1  and  0.5  were  supplied  to  the  porous 
sand  bed.  Electric  mass  flow  controllers  were  used 
to  control  the  mass  flow  rates  of  both  methane 
and  air.  The  mixture  velocity  was  defined  as 
U  =  m/(p0  A g),  where  m ,  p,  and  Ag  are  mass  flow 
rate,  density  at  room  temperature,  and  effective 
area  of  gas,  respectively.  The  effective  area  of 
gas  was  defined  as  Ag=As,  where  A  and  s  are 
the  cross-sectional  area  of  the  quartz  tube  and 
the  porosity  of  the  bed.  Additionally,  the  effective 
area  of  solid  was  As  =  A(l—s).  The  porosity  was 
defined  as  the  ratio  of  free  volume  to  the  total 
bed  volume  and  estimated  to  be  0.25  by  assuming 
a  homogeneous  porous  bed  consisting  of  spheres 
with  a  mean  diameter  of  0.56  mm. 

Five  thermocouples  (K-type,  500  pm,  and 
sheath  type)  were  arranged  along  the  centerline 
of  the  porous  sand  bed  with  50  mm  intervals  from 
the  downstream  side  of  combustor  as  shown  in 
Fig.  1.  To  determine  the  flame  location  in  the  por¬ 


Fig.  2.  Initial  temperature  profile  of  porous  bed  burner. 


ous  bed,  one  thermocouple  was  traversed  through 
the  porous  sand  bed  with  the  help  of  a  micrometer 
traverse  (minimum  resolution  of  0.05  mm  for  the 
flame  front  measurement  in  the  present  study) 
and  the  temperature  distribution  was  obtained. 
The  flame  was  assumed  to  be  located  at  a  point 
where  peak  temperature  was  observed.  Measure¬ 
ment  of  the  stabilized  flame  location  was  conduct¬ 
ed  under  the  steady  state  conditions.  The  system 
was  confirmed  to  be  in  steady  state  when  the  value 
of  temperature  fluctuation  remained  within  3  K 
during  1  h. 

3.  Experimental  results 

The  stabilized  flame  location  was  measured  by 
tracing  the  peak  temperature  position  by  travers¬ 
ing  a  thermocouple  in  the  axial  direction.  Varia¬ 
tion  of  the  stabilized  flame  location  with  mixture 
velocity  for  different  equivalence  ratios  (<f>  =  0.5 
and  0.1)  is  shown  in  Fig.  3.  It  was  difficult  to 
determine  the  flame  position  at  very  low  mixture 
velocities  due  to  very  low  heat  release  rate  and 
the  absence  of  an  explicit  temperature  peak. 
Therefore,  a  different  criterion  was  used  to  deter¬ 
mine  the  flame  position  at  low  mixture  velocities. 
The  flame  front  was  assumed  to  be  located  at  a 
position  where  the  difference  between  the  initial 
(without  mixture  flow)  and  final  (with  mixture 
flow)  temperature  profiles  of  the  porous  bed  was 
observed  to  be  maximum.  Figure  3  shows  that 
the  stabilized  flame  location  exhibits  S-shaped 
behavior  for  f  =  0.5.  This  S-shaped  behavior 
can  be  divided  into  three  regimes,  namely  low, 
moderate,  and  high  velocity  regimes.  In  the  low 
and  high  velocity  regimes,  the  flame  front  moves 
downstream  with  the  increase  of  the  mixture 
velocity.  Contrary  to  this  tendency  in  the  low 
and  high  velocity  regimes,  the  opposite  tendency 
of  flame  front  movement  in  the  upstream  direc¬ 
tion  is  observed  for  the  moderate  velocity  regime. 
This  is  due  to  strong  heat  recirculation  through 
the  porous  media  [17].  The  increase  of  the 


Fig.  3.  Experimental  result  of  variation  of  the  stabilized 
flame  location  with  mixture  velocity. 
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Fig.  4.  Direct  photographs  of  flame  front  at  0  —  0.5  (a)  U  —  0.24  m/s,  (b)  U  =  0.48  m/s,  and  (c)  U  —  0.56  m/s. 


incoming  enthalpy  with  the  increase  of  the  mix¬ 
ture  flow  rate  increases  the  effect  of  heat  recircula¬ 
tion  up  to  a  certain  value  of  the  mixture  flow  rate, 
which  results  in  a  flame  shifts  in  the  upstream 
direction.  However,  with  a  further  increase  in 
the  mixture  flow  rate,  the  residence  time  decreases 
and  this  becomes  dominant.  Since  the  amount  of 
the  incoming  enthalpy  is  very  small  for  0  =  0.1, 
such  S-shaped  behavior  was  not  observed  and 
the  flame  front  shifts  monotonically  in  the  down¬ 
stream  direction  with  an  increase  of  the  mixture 
velocity.  For  0  =  0.5,  oscillatory  flame  behavior 
was  observed  in  the  high  velocity  regime  and  for 
mixture  velocity  greater  than  0.48  m/s.  The  flame 
oscillated  between  two  axial  locations  which  were 
measured  by  tracing  the  position  of  the  peak  tem¬ 
perature.  The  approximate  variation  of  the  flame 
location  in  the  axial  direction  is  shown  in  Fig.  3 
with  open  square  symbols.  A  detailed  description 
of  this  oscillatory  behavior  will  be  presented  later. 

Figure  4  shows  flame  front  photographs  for 
the  moderate  and  high  velocity  regimes.  No  image 
for  the  low  velocity  regime  is  shown  because  lumi¬ 
nosity  from  the  reaction  zone  in  this  regime  was 
much  weaker  than  the  radiative  emission  from 
the  sand  bed.  Figure  4a  shows  the  flame  front  in 
the  moderate  velocity  regime  (0  =  0.5, 
U  =  0.24  m/s).  The  flame  front  is  nearly  flat  and 
it  is  stable  in  this  regime.  On  the  other  hand,  an 
inclined  flame  front  is  observed  in  the  high  veloc¬ 
ity  regime  (0  =  0.5,  U—  0.56  m/s)  as  shown  in 
Fig.  4c.  Between  these  two  regimes,  an  interesting 
feature  of  the  flame  front  is  observed  where  it 
breaks  down  into  two  parts:  top  and  bottom  as 
shown  in  Fig.  4b.  The  collapse  of  the  flame  front 
may  be  related  to  the  transition  of  the  combustion 
mode  from  a  stable  to  an  oscillatory  regime. 

Figure  5  shows  variation  of  the  oscillatory  time 
periods  with  the  mixture  velocity.  It  is  interesting 
to  note  that  the  oscillatory  time  periods  are  on  the 
order  of  an  hour  and  gradually  increase  with  mix¬ 
ture  velocity.  This  oscillatory  flame  behavior  was 
confirmed  to  be  independent  of  apparatus  related 
factors  since  small  fluctuations  in  heater  tempera¬ 
ture  (within  ±0.4%)  have  no  correlation  with  tem¬ 
perature  histories  in  the  porous  bed.  Although 
many  researchers  [19-21]  have  reported  the  exis¬ 
tence  of  oscillatory  flame  behavior  in  porous  beds, 
a  long  time  period  (on  the  order  of  an  hour)  of 


Fig.  5.  Experimental  results  of  variation  of  oscillation 
period  with  mixture  velocity. 


flame  oscillations  has  not  been  reported  in  the  lit¬ 
erature.  This  oscillatory  behavior  with  a  long  time 
period  is  probably  related  to  the  inclined  flame 
front  in  the  porous  sand  bed.  Since  the  value  of 
Gr/ Re1  based  on  the  pore  size  and  heater  temper¬ 
ature  as  an  ambient  of  the  porous  bed  is  very 
small  (<Cl)  in  the  present  system,  buoyancy  may 
not  play  a  significant  role  in  these  instabilities 
such  as  inclined  flame  front  and  oscillatory  behav¬ 
ior.  Thus,  we  suspect  that  the  cause  of  these  insta¬ 
bilities  may  be  related  to  the  diffusive-thermal 
instability,  including  large  radiative  heat  transfer. 
To  clarify  this,  further  studies  are  being  carried 
out  on  these  flame  oscillations. 


4.  Computation  model 

In  order  to  further  understand  the  experimen¬ 
tal  results,  one-dimensional  computational  studies 
were  carried  out  by  using  a  PREMIX  [22]-based 
code.  The  experimental  geometry  and  initial  tem¬ 
perature  profile  shown  in  Fig.  2  were  used  for  the 
present  calculations.  All  the  calculations  were 
conducted  for  methane-air  mixtures.  Energy 
equations  for  the  solid  and  gas  phases  are  as 
follows: 


■)  +—'52(>Y*v*c  P* 

/  cp  k=i 


dTg 

dx 
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+  —J2^khkWk+^hes{Tg 

CP  k=  1  CP 


Ts)  =  0 


(GasPhase) 


0 

dx 


+Ashgs(rs 


8=0 


(SolidPhase) 


Here, 

,  d Ts  ,  16a T\ 

q'  = kr  =  ~3P~ 


Q  =  axo(t*s  —  t\) 


The  effect  of  solid  conduction,  convective  heat 
transfer,  and  radiative  heat  transfer  in  the  porous 
sand  bed  was  considered.  Gas  phase  thermochem¬ 
ical  and  transport  properties  were  obtained  from 
the  Chemkin  [23]  and  Tranfit  [24]  packages.  The 
C-l  chemical  kinetics  mechanism  [25]  was  used. 
The  solid-to-solid  radiation,  hereafter  described 
as  solid  radiation,  was  approximated  using  opti¬ 
cally  thick  approximation  [26].  The  radiative 
properties  such  as  albedo,  emissivity,  and  extinc¬ 
tion  coefficient  were  assumed  to  be  unity.  The  heat 
transfer  coefficient  for  forced  convection  through 
the  packed  bed  [27]  was  calculated  by  using  the 
following  relation: 

Nu  =  .  fsDp,  ,  =  2A9(RePr)ip 
A(l-e)t 1/  v  ’ 


Here,  Dv,  2,  and  e  are  particle  diameter  of  spheres, 
thermal  conductivity  of  gas,  and  porosity,  respec¬ 
tively.  The  quantity  0  is  a  particle  shape  factor. 
Thermal  conductivity  of  solid  phase,  2S  was 
assumed  to  be  constant  [28].  Computations  for 
flame  behavior  such  as  the  flame  stabilization 
characteristic  were  conducted  to  interpret  the  ob¬ 
served  experimental  results  in  detail.  In  addition, 
the  effect  of  solid  radiation  on  the  stabilized  flame 
location  was  examined  by  comparing  the  two 
computational  models,  namely  the  case  with  solid 
radiation  (Model  I)  and  that  without  solid  radia¬ 
tion  (Model  II). 


5.  Computational  results  and  discussion 

5.7.  Stabilized  flame  location 

In  order  to  understand  the  various  combustion 
characteristics  of  gas  combustion  such  as  flame 
stabilization  and  its  structure  in  a  porous  sand 
bed,  computational  studies  were  carried  out. 
Figure  6  shows  the  computational  results  of 
variations  of  the  stabilized  flame  location  with 
mixture  velocity  at  0  =  0.1  and  0  =  0.5.  For 
purpose  of  comparison,  the  experimental  results 
are  also  presented  in  Fig.  6.  Computational  results 
without  solid  radiation  predicts  S-shaped 
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Fig.  6.  Computational  results  of  variation  of  flame 
location  with  mixture  velocity.  Here,  Model  I  includes 
solid  radiation  term  and  Model  II  is  without  solid 
radiation  term. 


behavior  for  the  stabilized  flame  location  at 
0  =  0.1,  while  those  with  solid  radiation  predict 
a  monotonic  dependence  on  the  stabilized  flame 
location,  which  is  similar  to  the  experimental 
results  for  0  =  0.1.  On  the  other  hand,  there  was 
a  limitation  point  for  the  steady  solution  in  the 
0  =  0.5  condition.  The  present  code  could  not 
predict  the  steady  solution  beyond  U—  0.372  m/s 
for  the  case  with  solid  radiation  and 
U—  0.167  m/s  for  that  without  solid  radiation. 
This  may  indicate  the  existence  of  unstable  solu¬ 
tions  such  as  oscillatory  behavior  beyond  these 
velocity  conditions.  Both  computational  results 
with  and  without  solid  radiation  show  a  trend 
similar  to  that  of  the  experimental  results.  It  is 
also  seen  that  the  prediction  results  with  solid 
radiation  (Model  I)  are  in  better  agreement  with 
the  experimental  results  than  the  predictions  from 
without  solid  radiation  (Model  II)  qualitatively. 

5.2.  Flame  structures 

The  flame  structures  in  the  moderate  and  high 
velocity  regimes  were  similar  to  a  conventional 
premix  flame.  The  typical  flame  structure  in  the 
high  velocity  regime  is  shown  in  Fig.  7a,  which 
shows  a  distinct  temperature  peak  along  with  a 
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Fig.  7.  Computational  result  of  the  flame  structure  (a) 
for  0  =  0.5,  U—  0.16  m/s,  (b)  for  0  =  0.5,  U  —  0.002 
m/s  with  solid  radiation. 

thin  reaction  zone.  The  methane  is  rapidly  con¬ 
sumed  near  the  reaction  zone  and  the  other  inter¬ 
mediate  chemical  species  such  as  CH,  CH3,  and 
CO  are  generated  within  the  thin  reaction  zone. 
However,  the  flame  structure  in  the  low  velocity 
regime  is  remarkably  different  from  that  in  the 
high  velocity  regime,  as  shown  in  Fig.  7b.  There¬ 
fore,  it  is  important  to  understand  the  flame  struc¬ 
ture  in  the  low  velocity  regime.  A  remarkable 
feature  of  the  flame  structure  in  the  low  velocity 
regime  is  the  absence  of  the  peak  temperature  in 
the  reaction  zone  and  existence  of  a  wide  reaction 
zone  as  indicated  by  the  width  of  CH3,  CH,  and 
other  species  profiles  in  Fig.  7b.  Since  no  peak 
temperature  in  the  reaction  zone  was  observed, 
the  location  of  the  maximum  CH  peak  in  the  cal¬ 
culation  was  considered  to  be  the  stabilized  flame 
location.  Despite  the  absence  of  a  temperature 
peak  in  the  reaction  zone,  combustion  reaction 
was  completed  as  shown  in  Fig.  7b.  CH4  is  grad¬ 
ually  consumed  throughout  the  reaction  zone. 
The  intermediate  species  such  as  CH3  and  CO 
are  generated  over  a  wider  zone,  which  is  approx¬ 
imately  100  mm.  Similar  combustion  phenomena 
with  a  wider  reaction  zone  are  observed  in  mild 
combustion  or  flameless  combustion  systems 
[29-31].  This  type  of  flameless  combustion 
behavior  in  the  low  velocity  regime  for  porous 
bed  combustion  has  not  been  reported  earlier. 


5.3.  Effect  of  solid  radiation  on  stabilized  flame 
location 

Figure  8  shows  the  variation  of  the  peak  tem¬ 
perature  in  the  reaction  zone.  The  peak  tempera¬ 
ture  increases  rapidly  up  to  a  mixture  velocity  of 
around  0.025  m/s,  and  beyond  that,  it  increases 
gently  with  the  increase  of  mixture  velocity.  Both 
computational  results  showed  the  same  trend.  On 
the  other  hand,  it  is  seen  that  the  difference  of 
peak  temperature  of  these  computation  results  is 
very  small  (~15  K  at  maximum). 

Figure  9  shows  a  comparison  of  the  computa¬ 
tional  flame  structure  at  0  =  0.5  and  U  —  0.16  m/s 
with  and  without  solid  radiation.  As  shown  in 
Fig.  9,  the  stabilized  flame  location  is  remarkably 
different  for  these  two  cases,  even  though  the 
results  were  obtained  at  the  same  equivalence 
ratio  and  velocity  condition.  Both  computations 
predicted  a  difference  of  peak  temperature  of  only 


Fig.  8.  Computational  results  of  variation  of  peak  flame 
temperature  with  mixture  velocity  for  0  =  0.5.  Model  I 
includes  solid  radiation  term  and  Model  II  is  without 
solid  radiation  term. 


Fig.  9.  Computational  result  of  comparison  of  flame 
structure  in  high  velocity  regime  for  0  =  0.5, 
U—  0.16  m/s.  Model  I  includes  solid  radiation  term 
and  Model  II  is  without  solid  radiation  term. 
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15  K  in  the  reaction  zone  as  shown  in  Fig.  8. 
However,  the  high  temperature  zone  generated 
by  chemical  reaction  for  the  case  of  without  solid 
radiation  is  much  broader  than  that  in  the  case 
with  solid  radiation  as  shown  in  Fig.  9.  In  addi¬ 
tion,  the  temperature  distribution  in  the  axial 
direction  for  the  case  with  solid  radiation  is  flatter 
than  the  temperature  distribution  for  the  case 
without  solid  radiation.  This  indicates  that  the 
heat  loss  from  the  high  temperature  zone  to  the 
lower  temperature  zone  upstream  and  down¬ 
stream  of  the  reaction  zone  is  much  higher  for 
the  case  with  solid  radiation  compared  with  that 
for  the  case  without  solid  radiation.  Therefore, 
the  total  heat  recirculation  rate  is  much  larger 
for  the  case  without  solid  radiation  when  com¬ 
pared  with  the  case  with  solid  radiation.  There¬ 
fore,  the  flame  front  is  stabilized  in  the  upstream 
region  in  the  case  without  solid  radiation. 

6.  Conclusions 

Combustion  characteristics  of  a  lean  methane- 
air  mixture  in  a  heated  porous  sand  bed  were 
investigated  both  experimentally  and  numerically. 
It  was  observed  that  the  stabilized  flame  location 
exhibits  S-shaped  behavior  for  0  =  0.5  with  mix¬ 
ture  velocity.  This  S-shaped  behavior  can  be  cate¬ 
gorized  into  three  regimes,  namely  high, 
moderate,  and  low  velocity  regimes.  For  the  mod¬ 
erate  velocity  regime,  a  stable  flat  flame  was 
observed.  Oscillatory  behavior  was  observed  in 
the  high  velocity  regime.  The  flame  front  oscillat¬ 
ed  with  a  characteristic  time  period  on  the  order 
of  lh  and  increased  with  the  increase  of  the  mix¬ 
ture  velocity.  An  inclined  flame  front  was 
observed  in  such  conditions.  A  simple  monoto¬ 
nous  dependence  of  flame  location  on  mixture 
velocity  was  observed  for  0  =  0.1.  The  effect  of 
solid  radiation  on  the  flame  behavior  was  exam¬ 
ined  through  computational  studies.  It  was  clari¬ 
fied  that  solid  radiation  plays  a  significant  role 
in  flame  stabilization.  From  the  examination  of 
the  flame  structure,  it  was  found  that  the  flame 
behavior  in  the  low  velocity  regime  was  similar 
to  that  in  the  flameless  combustion  mode. 
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Comments 

A.N.  Hay  hurst,  University  of  Cambridge ,  UK. 
Combustion  in  a  hot  bed  of  sand  is  inhibited  by 
the  recombination  of  radicals  on  the  exterior  of 
each  sand  particle  [1].  Because  of  this,  did  you 
vary  the  size  of  the  sand  particles  in  your  experi¬ 
ments  or  include  these  heterogeneous  effects  in  your 
modeling? 


Reference 

[1]  J.S.  Dennis,  A.N.  Hayhurst,  I.G.  Mackley,  Proc. 

Combust.  Inst.  19  (1982)  1205-1212. 

Reply.  We  did  not  vary  the  size  of  the  sand  particles 
in  the  present  study.  Kim  et  al.  [1]  showed  the  homoge¬ 
neous  chemical  reactions  overcome  the  effect  of  radical 
removal  by  the  surface  and  flame  becomes  resistant  to 
quenching  when  the  temperature  of  the  solid  wall  is 
above  873  K.  Since  mean  temperature  (1330  K)  of 
present  experiment  is  sufficiently  higher  than  the  solid 
temperature  (873  K),  we  did  not  consider  the  surface 
reaction  (radical  quenching  /  recombination  in  the 
present  work. 

We  did  not  include  the  heterogeneous  effects  into 
the  modeling  of  chemical  reaction  in  the  present 
computations  for  the  reasons  stated  above. 


Reference 

[1]  Kyu  Tae  Kim,  Dae  Hoon  Lee,  Sejin  Kwon,  Combust. 
Flame  146  (2006)  19-28. 


David  Smith,  University  of  Leeds,  UK.  Following  Al¬ 
lan  Hayhurst’s  comment,  as  an  alternative  to  changing 
particle  size  to  assess  the  significance  of  radical  loss  on 
sand  surfaces,  one  might  treat  the  particles  with  different 
wash  coatings  e.g.,  KC1  or  B203.  Classic  work  on  the  H2/ 
02  system  by  Baldwin  and  Walker  [1]  showed  that  this 
can  cause  very  different  surface  activity. 
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Reply.  The  effects  of  radical  recombination  and  rad¬ 
ical  quenching  are  not  significant  in  high  temperature 
(1330  K)  for  CH4/air  system  as  shown  in  Kim  et  al. 
([1]  in  reply  to  Hayhurst).  Therefore,  the  effect  of  radical 
loss  on  the  sand  surfaces  can  be  neglected  in  the  present 
study.  The  H2/02  system  might  be  active  at  lower  tem¬ 
perature,  since  the  sticking  coefficient  of  H2  is  much 
higher  than  that  of  CH4. 


